Improvement of risk estimate on wind turbine tower buckled by hurricane 
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Wind is one of the important reasonable resources. However, wind turbine towers are sure to be 
threatened by hurricanes. In this paper, method to estimate the number of wind turbine towers 
that would be buckled by hurricanes is discussed. Monte Carlo simulations show that our method 
is much better than the previous one. Since in our method, the probability density function of 
the buckling probability of a single turbine tower in a single hurricane is obtained accurately but 
not from one approximated expression. The result in this paper may be useful to the design and 
maintenance of wind farms. 

I. INTRODUCTION 

There arc rich offshore wind resources all over the world. Wind is the one with the largest installed-capacity growth 
from 2007 to 2012 among all renewable resources in America. U.S. wind power capacity increased from 8.7 GW 
in 2005 to 39.1 GW 2010 [1]. The National Renewable Energy Laboratory (NREL) estimates that offshore wind 
resources can be as high as four times the U.S. electricity generating capacity in 2010 [2]. 

Although the offshore wind resources are great, it is necessary to foresee the hurricane risks to offshore wind turbines. 
U.S. offshore resources are geographically distributed through the Atlantic, Pacific and Great Lake coasts. The most 
accessible shallow resources are located in the Atlantic and Gulf Coasts. Resources at depths shallower than 60 m in 
the Atlantic coast, from Georgia to Maine, are estimated to be 920 GW; the estimate for these resources in the Gulf 
coast is 460 GW [2]. Offshore wind turbines in these areas will be at risk from Atlantic hurricanes. Between 1949 
and 2006, 93 hurricanes struck the U.S. mainland according to the HURDAT (Hurricane Database) database of the 
National Hurricane Center [.')]. Hurricane risks are quite variable, both geographically and temporally. Pielke et al. 
note pronounced differences in the total hurricane damages (normalized to 2005) occurring each decade [4]. Hurricanes 
and other hazards can cause widespread electric power outages, which, in turn, can affect business operations, heating, 
financial transactions, security systems, water distribution, traffic signaling, and countless other aspects of daily life. 
And hurricane hazards cause the most damage to power systems in the Eastern U.S ["i]. Interdecadal major hurricane 
fluctuations occur in both landfall locations and overall activity [G-9] . Most of the deadliest and costliest Atlantic 
tropical cyclones are major hurricanes [C]. Major hurricanes account for just over 20% of the tropical storms and 
hurricanes that strike the United States but cause more than 80 percent of the damage [10]. 

Hurricanes Katrina and Rita hit the center of the American petrochemical industry, shutting down eight refineries, 
hundreds of oil-drilling and production platforms, and numerous other industrial facilities. Furthermore, they triggered 
numerous hazardous-materials (hazniat) releases from industrial facilities and storage terminals onshore, as well as 
from oil and gas production facilities offshore in the Gulf of Mexico (GoM) [U]. Hurricane Ivan had caused much 
concern among industrialists, operators, and government officials on the performance of the offshore oil and gas drilling 
and production activities and infrastructure in the GoM during major hurricanes [12]. Hurricanes in the GoM, which 



can reach wind speeds of 240 km/h accompanied by waves of over 25 m, pose serious challenges for the design and 
operation of offshore facilities in these harsh environments [1-3]. 

The 2005 Atlantic hurricane season was the most active hurricane season on record with 28 named storms (previous 
record was 21 named storms set in 1933), 15 of which reached hurricane status (previous record was 12 set in 1969) 
[14]. The IPCC suggests in 1990 that: there are some evidences from model simulations and empirical considerations 
that the frequency per year, intensity and area of occurrence of tropical disturbances may increase in a doubled carbon 
dioxide world, though it is not compelling [15]. Recent sizable hurricane losses have raised questions about the causes. 
Some have claimed that storm frequencies and/or intensities have increased [l(i, 17], but other studies indicate no 
long-term trend in hurricane activity during the 20th century [In]. Others see the increases as a result of climate 
change resulting from global warming [19, 20]. Others consider the higher hurricane losses a result of societal changes 
leading to greater vulnerability in hurricane prone areas [21]. 

The average annual insured losses from hurricanes are 2.6 billion dollars for 1949-2006, the highest storm-related 
loss in the nation. The ranks of the average annual insured property losses of all the extreme weather conditions in the 
U.S. (2006 dollars in billions) are: (1) hurricanes (2.6 billion dollars), (2) floods (2.2 billion dollars), (3) thunderstorms 
(1.6 billion dollars), (4) tornadoes (1.0 billion dollars), (5) hail (0.9 billion dollars), (6) snowstorms (0.5 billion dollars), 
(7) freezing rainstorms (0.2 billion dollars), and (8) wind storms (0.2 billion dollars) [22]. On the other hand, world 
energy problems are becoming more and more serious. It is imminent to make use of rich wind resources. While the 
development of onshore wind turbines has gone on wheels, there are 20 offshore wind projects in the planning process 
[2]. 

So it is important to analyze the hurricane risks to offshore wind turbines. The design, maintenance and assessment 
of offshore structures must meet the requirements laid down in the Code of Federal Regulations, Title 30, Part 250 [2:-]]. 
And the design requirements for offshore facilities follow the recommendations of the American Petroleum Institute 
which calls for platforms and floating permanent systems to be designed to withstand a full-population hurricane [24] 
with a return period of 100 years [25, 20]. These 100-year criteria correspond to a wind speed of about 150 km/h in 
1-h average winds or about 180 km/h in sustained 1-min winds and a maximum wave height of 22m [12, 27, 28]. 

In [29], a series of models used to describe the characters of hurricanes and turbine towers are established, based on 
which the risks on offshore wind turbines from hurricanes are discussed, and a model used to estimate them in four 
representative locations (Galveston County, TX; Dare County, NC; Atlantic County, NJ; and Dukes County, MA) in 
the Atlantic and Gulf Coastal waters of the United States has been built. Although the models in [2!)] can already 
give a good and reasonable risk estimate of hurricanes, one approximated expression used to calculate the probability 
density function of turbine buckling probability is able to be improved. In this study, we will give out an accurate 
expression of it. Some discussions and results given in [29] are based on Monte Carlo simulations. One of them has its 
mathematical expression which will be discussed later. Thanks to this work, Monte Carlo simulations can be avoided. 
Perhaps it will be helpful for further theoretical analyses. At last, we will put forward an idea of estimating hurricane 
risks which depends on the expected survive time (EST). 

The organization of this paper is as follows. The basic models and the improved methods will be briefly described 
in the next section, and then in Sec. Ill the analyses on Monte Carlo simulations will be presented. In sec. IV, 
a short discussion about expected survive time (EST) will be given. The results will be given in Sec. V, and then 
concluding remarks in the final section. 



II. MODELS OF HURRICANES AND TURBINES 

A. The models established in [29] 

The previous model presented by Stephen et al. in [29] can be summarized as follows. 

(1) Hurricane occurrence. 

Hurricane occurrence is modeled as a Poisson process with rate A obtained by fitting to historical hurricane data. 
The probability that H, the number of hurricanes that occur in T years, equals a particular value h is given by: 

PriH = h) = ^e-^^. (1) 

(2) The maximum 10-min sustained wind speed of each hurricane at 10-m height. 

It is supposed that there is a maximum 10-min sustained wind speed during each hurricane. The wind speed at 
90-m height (hub- height of turbine towers) decides the probability of a single wind turbine tower buckling [30] . The 
maximum 10-min sustained wind speed of each hurricane at 10-m height is modeled as Generalized Extreme Value 
(GEV) distribution with a location parameter ^, a scale parameter a, and a shape parameter ^ fitted to historical 
hurricane data. The probability density function for W, the maximum sustained wind speed, evaluated at particular 
value w is given by: 




/H^H==-cxp - 1 + e — ^ i + e^^-^ • (2) 



(3) The buckling probability of a single wind turbine tower. 

The probability that a single wind turbine tower is buckled by a maximum 10-min sustained hub-height wind speed 
u is modeled using a log- logistic function with a scale parameter a and a shape parameter 13. These parameters 
are fitted to probabilities of turbine tower buckling calculated by comparing the results of simulations of the 5-MW 
offshore wind turbine designed by the NREL to the stochastic resistance to buckling proposed by S0rensen, et al. 
[31, 32]. The buckling probability of a single wind turbine B can be given by the following log-logistic function: 

B{u) = ^"/ (3) 

In order to distinguish the sign of buckling probability from the differential sign, we use B to express the random 
variable and b to express its exact value instead of D and d used by Stephen et al. originally. 

(4) Fitting to Beta distribution. 

The exact hub-height (90-m) wind speed u is related to the exact value w of 10-m height random wind speed W . 
This wind speed is scaled from 90-m height to 10-m height assuming power-law wind shear with an exponent of 0.077 
[33]. It can be expressed as: 

Iniulw) , , 

077 = — — — — f4) 

ln(90/10) ^^ 

A Beta distribution is fitted to Monte Carlo simulations of the convolution of B and W . The procedure used in 
[29] is as follows: 
1. Simulate a large number of wind speeds W with a Generalized Extreme Value (GEV) distribution in Eq. (2). 



2. Calculate the probability of turbine tower buckling B for each wind speed from step. 1 using the log-logistic 
damage function in Eq. (3) and the wind-shear function in Eq. (4). 

3. Calculate the empirical cumulative distribution function (CDF) for the buckling probabilities in step. 2. For each 
probability of buckling, calculate the probability that value occurs. The result is an x-y graph with Probability of 
turbine tower buckling on the x-axis and Probability of occurrence on the y-axis. 

4. Use nonlinear curve fitting to fit the CDF of a Beta distribution to the empirical CDF in step. 3. Use starting 
values of as ~ 0.02 and 13b = 0.2. 

The probability density function (PDF) of beta distribution with parameters as and I3b is given by: 

f{x;aB,l3B)^-l — — — -• (5) 

(5) The probability of buckling turbine number. 

Through Eq. (5), it is easy to show that the number of turbines buckled by a single hurricane in a wind farm with 
n turbines can be modeled by a beta-binomial distribution with parameters as and /3b- The probability that X, the 
number of turbine towers that buckle out of n total, equals a particular value x is given by: 

Pr{X^x)^[ — -— . (6) 

\xj B(aB,PB) 

(6) The number of buckling turbines in time period T. 

The cumulative distribution of the buckling turbine number in T years without replacement, Yno rep, is modeled as 
a modified phase- type distribution [34, 35]: 

Pr{Yno rop <y\T<t)^ gexp(TT(y, n))e. (7) 

where g, e are vectors, and T is one matrix (see [29] for detailed explanation for them). 

B. Our improved method 

Step. (l)-(3) are the methods to establish basic mathematical models of hurricanes and turbine towers. Step. 
(4)- (6) are the methods to solve the models. 

Monte Carlo simulations are used to estimate the CDF of the buckling probability of a single turbine tower in a 
single hurricane through the procedure presented in step. (4). In fact, the buckling probability for a single turbine 
in a hurricane _B is a function of maximum 10-min sustained wind speed u at hub- height, while u and w fit Eq. (4), 
so its PDF can be calculated from the PDF of W which is given by Eq. (2). Here we need a theorem: 
Theorem if F = f{X) is a monotonic function, X = /"^(F) £ C^{R), and the PDF of X is p{x), then the PDF of 
y, q{y) is given by: 

q{y)^p{f-Hy))\^^^l (8) 

■ 

By taking use of the above theorem, we can get the PDF of B noted as p{b) : 

p& =-exp-l + e Ol + C «-^- -T-7, 9 



where it is supposed that u — sw. s ~ exp(0.0771og(9)) can be got from Eq. (4). 

In this paper, Eq. (9) will be used to estimate the number of buckling turbines in time period T directly. Let A{h) 
be one (n + 1) x (71 + 1) matrix with its elements given by: 




A{h){i,j) = { ' '-^ ' ' ■ - , (10) 



where h is the buckling probability of a single wind turbine after a single hurricane and n is the number of turbines in 
one farm. We also define one state vector f = (/i, /2, fs,- ■ ■ , fn+i) with fi being the probability that i — 1 turbines 
are buckled out of n total after a single hurricane, f = g = (1, 0, 0, ■ • • , 0) at first, since none turbine is buckled. It 
is obvious that the state vector f ~ sY[i=i -^i^j) a-fter k hurricanes, where bj is a sample of b in the jth hurricane. 
According to Eq. (1). the probability that k hurricanes happen in T years is ^ ^,' e^^^. The expected value E(i) of 
f can be obtained as follows: 

E{i) = E(e-'^ + f^sllAib,)i^e-A , (11) 



hi 

fc=i i=i 

where bj (j' = 1, 2, • • • ) are independent of each other according to the assumption in [29] that the maximum 10-min 
sustained wind speed of each hurricane is independent of each other. So £'(ni=i ^(^j)) = Ylj=i ^(^i^j))- Notice 
that E{A{bi)) = E{A{bj)) = E{A{b)), since bj (j = 1, 2, • • • ) are samples of b. Eq. (11) is equivalent to: 

E{i) = gexp(Ar£;(A(6)))e-^^ ^ gcxp{XT{E{A{b)) - /)). (12) 

Compare it with the expression gexp(TT(j/,n)) of f in Eq. (7), they are same except replacing the matrix T in Eq. 
(7) by X{E{A{b)) — /). It can be expressed as follows: 

where B{) is the beta function in Eq. 6. And we have: 

E (('' T^ j ^]&J-*(1 - b)"-^+A = r ('' T^l ^^-'(1 - br-^+'pib)db. (14) 

ai and 02 are the bounds of b which are involved by the GEV distribution of W. As mentioned before, p{b) is directly 
used in this paper. Since the essence of this method is given out by Eq. (13), we can use it to calculate l^rcp which 
is a estimate of the buckling turbine number in T years with replacement after each hurricane (see [29] for detailed 
explanation for it). 

It will be seen in See. V that this method gives a closer estimate to the Monte Carlo simulations in both cases with 
or without replacement. 



III. THE ANALYSIS ON MONTE CARLO SIMULATIONS 

In [29], the loss caused by the Category (a classification standard of hurricanes) 1 to 3 hurricanes is estimated by 
doing Monte Carlo simulations and excluding the simulations with Category 4 to 5 hurricanes happening manually. 



which is equivalent to replace the PDF of B by the condition PDF of B ioi b < fosmax- ^Smax is the maximum buckling 
probability caused by a Category 3 hurricane of a single turbine tower. By Eq. (3) and Eq. (4) one can get: 



JSmax o I 

1+ sm^^Y 



(15) 



where wsmax is the largest value of maximum 10-min sustained wind speed at 10-m height of Category 3 hurricanes. 
u — sw and s = cxp(0.0771og(9)) give the meaning of s. With the definition of b^n-iax, we have the condition PDF of 
B ior b < 63max as: 

It will be seen in Sec. V that this method fits the Monte Carlo simulations very well. 



IV. EXPECTED SURVIVE TIME (EST) 

In this section, one new way will be given to analyze the risks of a wind farm suffered from hurricanes. It is supposed 
that the buckling of each turbine is independent of each other. If the character of an individual turbine buckling is 
confirmed, it can describe the character of the whole wind farm in some aspects. The probability q of one turbine 
not buckling in io years is given by: 

k=ij=i 



which is equivalent to 



q = e-^*«^W. (18) 



Here bj means the probability of this individual turbine buckling in jth hurricane. Naturally the probability of one 
turbine buckling in to years, noted as 1 — q is given by: 

l-q=l-e-^*°^^''\ (19) 

Actually one turbine buckling in to years means the exact value t of its survive time T is less than to. So 1 — q is the 
CDF of T. The PDF of T, noted as p{t) can be expressed as: 

p{t) = XE{b)er^'^^''\ (20) 

We have the expected survive time E(T) as: 

EiT)^ftpm^^^. (21) 

The expected survive time E{T) means the mean value of survive time T. It is determined by £'(6) and A only. E{b) 
and A are parameters which reflect the surrounding model of the wind farm. The worse is the surrounding, the less 
is E{T). So E{T) can be used to compare wind farms with each other or decide whether a project should be applied 
to a wind farm. We have E{b) as: 

E{b) = / bp{b)db. (22) 



It is easy to show that: 

l-g=l-e~'^ (23) 

by using Eq. (19) and Eq. (21). Wc can get the probabihty of an individual turbine buckling in T years only with 
EST known. The expected survive number (ESN) E{1) can be calculated by: 

E{l)=n{l~q)^n{l-e~^)), (24) 

where n is the total turbine number. It will be seen in See. V that ESN got by Eq. (24) is close to the results of 
Monte Carlo simulations and the state vector calculated by the method introduced in Sec. II. 

V. RESULTS 

In our calculating, the same data and parameters in [29] are used (for details, see the captions of Figs. 1-3). Firstly, 
our model is tested in calculating the CDF of buckling number without replacement in the wind farm of Galveston 
County, TX. Suppose that turbines are pointed into wind (Active Yawing). Test period T and total turbine number 
n are set to be 20 and 50 respectively. So are the following tests. We always use full line in red for the new method, 
dotted line in blue for Monte Carlo simulations and chain line in black for the original method in [2!)]. The result of 
new method is almost same to Monte Carlo simulations according to Fig. (1), which shows the accuracy of the new 
method is better. 

Secondly, the new model is tested when buckling turbines are replaced after each hurricane in the same wind farm 
(Galveston County, TX). But this time, turbines are pointed perpendicular to wind (Not Yawing). Again the new 
method is closer to Monte Carlo simulations according to Fig. (2). 

In these two tests listed above, the new method gives more accurate results of the losses, which may be helpful for 
some further estimates and analyses on risks suffered from hurricanes. 

Then we test the condition PDF model mentioned in Sec. Ill in the wind farm of Dare County, NC. Suppose that 
turbines are pointed perpendicular to wind (Not Yawing). The boundary value of 10-m wind speed between Category 
3 and Category 4 hurricanes is not given out in [:il)]. So we choose it to be 113 knots. One can find in Fig. (3) that 
the condition PDF of B can give out a good estimate of the CDF of buckling number under a given condition, since 
its curve fits Monte Carlo simulations accurately. 

Finally, the expected survive number (ESN) is calculated in three different ways. They arc the state vector given 
by Eq. (12), Monte Carlo simulations based on Eq. (1), Eq. (2) and Eq. (3), the expected buckling probability given 
by Eq. (19) respectively. The wind farm is chosen in TX and turbines are active yawing. These three results are 
extremely close to each other, which means that the buckling probability of a single turbine may refiect the risks of 
a wind farm suffered from hurricanes there reasonably. 

As mentioned in See. Ill, the mean of state vector in T years without Category 4 to 5 hurricanes happening can 
be described easily by using a condition PDF oi B. It is different from the situation that Category 4 to 5 hurricanes 
happen between Category 1 to 3 hurricanes while only the damage caused by Category 1 to 3 hurricanes is considered. 
It is easy to analyze it by Monte Carlo simulations. We have tried to give out the mathematical expression of the state 
vector caused by Category 1 to 3 hurricanes only when Category 1 to 3 hurricanes and Category 4 to 5 hurricanes 



happen alternately. However, it is difficult to handle this because of the correlations between hurricanes. This work 
should be meaningful since we can ensure the percentage of buckling turbines of each hurricane Category by the 
parameters of the surrounding in a wind farm [such as A in Eq. (1)] and the turbine [sucn as a and /3 in Eq. (3)] 
directly but not by Monte Carlo simulations. 

In Sec. IV, the buckling probability of a single turbine in T years is proposed. We notice that the buckling of each 
turbine in each hurricane is independent of each other. So we try to get the state vector of turbine buckling number 
after T years by using Eq. (23). However, the turbines in a same wind farm are in fact not independent, because 
they are sure to suffer from same hurricanes during T years. It is not same to the case that each turbine is located 
in a different wind farm with same surrounding. So Eq. (23) can not take place of the state vector calculated by Eq. 
(12) completely. 

VI. CONCLUDING AND REMARKS 

In this paper, the characters of hurricanes and wind turbines are discussed by similar methods as established by 
Rose et. al in [2!)]. In which the risk of wind turbine towers suffered from hurricanes is estimated by cumulative 
distribution function (CDF) of the number of buckling turbine towers, and the hurricane risks in a wind farm are 
analyzed by the expected survived time (EST) of a single turbine. Monte Carlo simulations show that our results are 
accurate enough. The study in this paper is helpful to understand the effects of hurricanes on a wind farm, and may 
also be valuable to the design and maintenance of wind farms. 
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FIG. 1: The cumulative distribution function (CDF) of buckling number of 50 turbines totally in 20 years. Buckling turbines 
aren't replaced. a;-axis is the buckling number of wind turbine towers after 20 years. The value of y-axis at x is the probability 
that the buckling number of wind turbine towers is less than x. Turbines are actively yawing in Galveston County, TX. The 
CDF is obtained by Eq. (9) with parameters a = 174, /3 = 19.3, A = 0.19, p = 78.7, a = 12.1, ^ = 0.251. 
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FIG. 2: The cumulative distribution function (CDF) of buckling number of 50 turbines totally in 20 years. Buckling turbines 
are replaced after each hurricane, s-axis is the buckling number of wind turbine towers after 20 years. The value of y-axis 
at X is the probability that the buckling number of wind turbine towers is less than x. Turbines are not yawing in Galveston 
County, TX. The CDF is obtained by Eq. (9) with parameters a = 140, (3 = 18.6, A = 0.19, ^l = 78.7, a = 12.1, ^ = 0.251. 
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FIG. 3: The cumulative distribution function (CDF) of buckling number of 50 turbines totally in 20 years without Category 4 
to 5 hurricanes. Buckling turbines aren't replaced, x-axis is the buckling number of wind turbine towers after 20 years caused 
by Category 1 to 3 hurricanes. The value of y-axis at x is the probability that the buckling number of wind turbine towers is 
less than x. Turbines are not yawing in Dare County, NC. The CDF is obtained by Eq. (16), Eq. (9) with parameters a = 140, 
/3 = 18.6, A = 0.21, n = 77.6, a = 11.9, ,^ = -0.0366 and Eq. (15) with parameters w = 113. 



Method 



State Vector 
Monte Carlo 
Expected buckling probability 



ESN 



5.8884 
5.8412 
5.8885 



TABLE I: The expected survive number (ESN) of 50 turbines totally in 20 years is calculated in three different ways. Turbines 
are actively yawing in Galveston County, TX. The parameter values used in these calculations are a — 174, /3 = 19.3, A = 0.19, 
^ = 78.7, a = 12.1, ^ = 0.251, see Eq. (2), (3) and (9). 



